MEMOR: A database of archeological human remains collections from Flanders, Belgium.

The aim of this article is to describe a newly created open access database of archeological human remains collections from Flanders, Belgium. The MEMOR database (www.memor.be) was created to provide an overview of the current practices of loans, reburial, and the research potential of human skeletons from archeological sites currently stored in Flanders. In addition, the project aimed to provide a legal and ethical framework for the handling of human remains and was created around stakeholder involvement from anthropologists, geneticists, contract archeologists, the local, regional and national government agencies, local and national government, universities, and representatives of the major religions. The project has resulted in the creation of a rich database with many collections available for study. The database was created using the open-source Arches data management platform that is freely available for organizations worldwide to configure in accordance with their individual needs and without restrictions on its use. Each collection is linked to information about the excavation and the site the remains originate from, its size and time period. In addition, a research potential tab reveals whether any analyses were performed, and whether excavation notes are available with the assemblage. The database currently contains 742 collections, ranging in size from 1 to over 1000 individuals. New collections will continue to be added when new assemblages are excavated and studied. The database can also be expanded to include human remains collections from other regions and other material categories, such as archaeozoological collections.

Genomic analyses of hair from Ludwig van Beethoven.

Ludwig van Beethoven (1770-1827) remains among the most influential and popular classical music composers. Health problems significantly impacted his career as a composer and pianist, including progressive hearing loss, recurring gastrointestinal complaints, and liver disease. In 1802, Beethoven requested that following his death, his disease be described and made public. Medical biographers have since proposed numerous hypotheses, including many substantially heritable conditions. Here we attempt a genomic analysis of Beethoven in order to elucidate potential underlying genetic and infectious causes of his illnesses. We incorporated improvements in ancient DNA methods into existing protocols for ancient hair samples, enabling the sequencing of high-coverage genomes from small quantities of historical hair. We analyzed eight independently sourced locks of hair attributed to Beethoven, five of which originated from a single European male. We deemed these matching samples to be almost certainly authentic and sequenced Beethoven's genome to 24-fold genomic coverage. Although we could not identify a genetic explanation for Beethoven's hearing disorder or gastrointestinal problems, we found that Beethoven had a genetic predisposition for liver disease. Metagenomic analyses revealed furthermore that Beethoven had a hepatitis B infection during at least the months prior to his death. Together with the genetic predisposition and his broadly accepted alcohol consumption, these present plausible explanations for Beethoven's severe liver disease, which culminated in his death. Unexpectedly, an analysis of Y chromosomes sequenced from five living members of the Van Beethoven patrilineage revealed the occurrence of an extra-pair paternity event in Ludwig van Beethoven's patrilineal ancestry.

Large-scale pedigree analysis highlights rapidly mutating Y-chromosomal short tandem repeats for differentiating patrilineal relatives and predicting their degrees of consanguinity.

Rapidly mutating Y-chromosomal short tandem repeats (RM Y-STRs) were suggested for differentiating patrilineally related men as relevant in forensic genetics, anthropological genetics, and genetic genealogy. Empirical data are available for closely related males, while differentiation rates for more distant relatives are scarce. Available RM Y-STR mutation rate estimates are typically based on father-son pair data, while pedigree-based studies for efficient analysis requiring less samples are rare. Here, we present a large-scale pedigree analysis in 9379 pairs of men separated by 1-34 meioses on 30 Y-STRs with increased mutation rates including all known RM Y-STRs (RMplex). For comparison, part of the samples were genotyped at 25 standard Y-STRs mostly with moderate mutation rates (Yfiler Plus). For 43 of the 49 Y-STRs analyzed, pedigree-based mutation rates were similar to previous father-son based estimates, while for six markers significant differences were observed. Male relative differentiation rates from the 30 RMplex Y-STRs were 43%, 84%, 96%, 99%, and 100% for relatives separated by one, four, six, nine, and twelve meioses, respectively, which largely exceeded rates obtained by 25 standard Y-STRs. Machine learning based models for predicting the degree of patrilineal consanguinity yielded accurate and reasonably precise predictions when using RM Y-STRs. Fully matching haplotypes resulted in a 95% confidence interval of 1-6 meioses with RMplex compared to 1-25 with Yfiler Plus. Our comprehensive pedigree study demonstrates the value of RM Y-STRs for differentiating male relatives of various types, in many cases achieving individual identification, thereby overcoming the largest limitation of forensic Y-chromosome analysis.

From collected stamps to hair locks: ethical and legal implications of testing DNA found on privately owned family artifacts.

Biological samples containing DNA that is attributed to deceased relatives, can now undergo genetic testing at a reasonable cost due to revolutionary improvements in sampling, sequencing, and analytical techniques. This artifact DNA testing, or 'artDNA', includes genetic analysis of hair locks, stamps, envelopes with saliva traces or teeth. ArtDNA can reveal valuable information about a deceased relative or one's genetic background, but it also presents novel ethical dilemmas and legal uncertainties for genetic researchers and commercial testing services. In this paper, we provide an analysis of some of the unique ethical and legal risks of such testing and provide needed recommendations for practitioners of private family artDNA testing. ArtDNA testing generates ethical and legal risks regarding the privacy and autonomy of deceased individuals, the rights of living relatives over their ancestor's genetic information, and the rights of living persons to control their own genetic information. To mitigate these risks, practitioners can conduct certain preliminary testing to ascertain the identity of a DNA donor and estimate the time that has elapsed postmortem. Generally, the ethical and legal concerns will be higher when a shorter period has passed between the death of the DNA donor and the time of artifact DNA testing. Regardless, all artDNA testing present some risks, and practitioners should exercise professional judgement as necessary.

Y-chromosome target enrichment reveals rapid expansion of haplogroup R1b-DF27 in Iberia during the Bronze Age transition.

The Y chromosome can yield a unique perspective into the study of human demographic history. However, due to the repetitive nature of part of its sequence, only a small set of regions are suitable for variant calling and discovery from short-read sequencing data. These regions combined represent 8.9 Mbp or 0.14% of a diploid human genome. Consequently, investing in whole-genome sequencing to resolve Y-chromosome questions is poorly efficient. Here we use, as an alternative, target enrichment technology to greatly increase sequencing effectiveness, validating and applying the technique to 181 males, for 162 of whom we obtained a positive result. Additionally, 75 samples sequenced for the whole genome were also included, for a total sample size of 237. These samples were chosen for their Y chromosome haplogroup: R1b-DF27. In the context of European populations, and particularly in Iberia, this haplogroup stands out for its high frequency and its demographic history. Current evidence indicates that the diffusion of this haplogroup is related to the population movements that mark the cultural Bronze Age transition, making it remarkably interesting for population geneticists. The results of this study show the effects of the rapid radiation of the haplogroup in Spain, as even with the higher discriminating power of whole sequences, most haplotypes still fall within the R1b-DF27* paragroup rather than in the main derived branches. However, we were able to refine the ISOGG 2019-2020 phylogeny, and its two main subbranches, namely L176.2 and Z272, which present geographical differentiation between the Atlantic and Mediterranean coasts of Iberia.

3D facial phenotyping by biometric sibling matching used in contemporary genomic methodologies.

The analysis of contemporary genomic data typically operates on one-dimensional phenotypic measurements (e.g. standing height). Here we report on a data-driven, family-informed strategy to facial phenotyping that searches for biologically relevant traits and reduces multivariate 3D facial shape variability into amendable univariate measurements, while preserving its structurally complex nature. We performed a biometric identification of siblings in a sample of 424 children, defining 1,048 sib-shared facial traits. Subsequent quantification and analyses in an independent European cohort (n = 8,246) demonstrated significant heritability for a subset of traits (0.17-0.53) and highlighted 218 genome-wide significant loci (38 also study-wide) associated with facial variation shared by siblings. These loci showed preferential enrichment for active chromatin marks in cranial neural crest cells and embryonic craniofacial tissues and several regions harbor putative craniofacial genes, thereby enhancing our knowledge on the genetic architecture of normal-range facial variation.

Subdividing Y-chromosome haplogroup R1a1 reveals Norse Viking dispersal lineages in Britain.

The influence of Viking-Age migrants to the British Isles is obvious in archaeological and place-names evidence, but their demographic impact has been unclear. Autosomal genetic analyses support Norse Viking contributions to parts of Britain, but show no signal corresponding to the Danelaw, the region under Scandinavian administrative control from the ninth to eleventh centuries. Y-chromosome haplogroup R1a1 has been considered as a possible marker for Viking migrations because of its high frequency in peninsular Scandinavia (Norway and Sweden). Here we select ten Y-SNPs to discriminate informatively among hg R1a1 sub-haplogroups in Europe, analyse these in 619 hg R1a1 Y chromosomes including 163 from the British Isles, and also type 23 short-tandem repeats (Y-STRs) to assess internal diversity. We find three specifically Western-European sub-haplogroups, two of which predominate in Norway and Sweden, and are also found in Britain; star-like features in the STR networks of these lineages indicate histories of expansion. We ask whether geographical distributions of hg R1a1 overall, and of the two sub-lineages in particular, correlate with regions of Scandinavian influence within Britain. Neither shows any frequency difference between regions that have higher (≥10%) or lower autosomal contributions from Norway and Sweden, but both are significantly overrepresented in the region corresponding to the Danelaw. These differences between autosomal and Y-chromosomal histories suggest either male-specific contribution, or the influence of patrilocality. Comparison of modern DNA with recently available ancient DNA data supports the interpretation that two sub-lineages of hg R1a1 spread with the Vikings from peninsular Scandinavia.

Identification and characterization of novel rapidly mutating Y-chromosomal short tandem repeat markers.

Short tandem repeat polymorphisms on the male-specific part of the human Y-chromosome (Y-STRs) are valuable tools in many areas of human genetics. Although their paternal inheritance and moderate mutation rate (~10-3 mutations per marker per meiosis) allow detecting paternal relationships, they typically fail to separate male relatives. Previously, we identified 13 Y-STR markers with untypically high mutation rates (>10-2 ), termed rapidly mutating (RM) Y-STRs, and showed that they improved male relative differentiation over standard Y-STRs. By applying a newly developed in silico search approach to the Y-chromosome reference sequence, we identified 27 novel RM Y-STR candidates. Genotyping them in 1,616 DNA-confirmed father-son pairs for mutation rate estimation empirically highlighted 12 novel RM Y-STRs. Their capacity to differentiate males related by 1, 2, and 3 meioses was 27%, 47%, and 61%, respectively, while for all 25 currently known RM Y-STRs, it was 44%, 69%, and 83%. Of the 647 Y-STR mutations observed in total, almost all were single repeat changes, repeat gains, and losses were well balanced; allele length and fathers' age were positively correlated with mutation rate. We expect these new RM Y-STRs, together with the previously known ones, to significantly improving male relative differentiation in future human genetic applications.

Correction to: The Dutch Y-chromosomal landscape.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

The Dutch Y-chromosomal landscape.

Previous studies indicated existing, albeit limited, genetic-geographic population substructure in the Dutch population based on genome-wide data and a lack of this for mitochondrial SNP based data. Despite the aforementioned studies, Y-chromosomal SNP data from the Netherlands remain scarce and do not cover the territory of the Netherlands well enough to allow a reliable investigation of genetic-geographic population substructure. Here we provide the first substantial dataset of detailed spatial Y-chromosomal haplogroup information in 2085 males collected across the Netherlands and supplemented with previously published data from northern Belgium. We found Y-chromosomal evidence for genetic-geographic population substructure, and several Y-haplogroups demonstrating significant clinal frequency distributions in different directions. By means of prediction surface maps we could visualize (complex) distribution patterns of individual Y-haplogroups in detail. These results highlight the value of a micro-geographic approach and are of great use for forensic and epidemiological investigations and our understanding of the Dutch population history. Moreover, the previously noted absence of genetic-geographic population substructure in the Netherlands based on mitochondrial DNA in contrast to our Y-chromosome results, hints at different population histories for women and men in the Netherlands.

Ysurnames? The patrilineal Y-chromosome and surname correlation for DNA kinship research.

The Y-chromosome is a widely studied and useful small part of the genome providing different applications for interdisciplinary research. In many (Western) societies, the Y-chromosome and surnames are paternally co-inherited, suggesting a corresponding Y-haplotype for every namesake. While it has already been observed that this correlation may be disrupted by a false-paternity event, adoption, anonymous sperm donor or the co-founding of surnames, extensive information on the strength of the surname match frequency (SMF) with the Y-chromosome remains rather unknown. For the first time in Belgium and the Netherlands, we were able to study this correlation using 2,401 males genotyped for 46 Y-STRs and 183 Y-SNPs. The SMF was observed to be dependent on the number of Y-STRs analyzed, their mutation rates and the number of Y-STR differences allowed for a kinship. For a perfect match, the Yfiler® Plus and our in-house YForGen kit gave a similar high SMF of 98%, but for non-perfect matches, the latter could overall be identified as the best kit. The SMF generally increased due to less mismatches when encountering [1] deep Y-subhaplogroups, [2] less frequently occurring surnames, and [3] small geographical distances between relatives. This novel information enabled the design of a surname prediction model based on genetic and geographical distances of a kinship. The prediction model has an area under the curve (AUC) of 0.9 and is therefore useable for DNA kinship priority listing in estimation applications like forensic familial searching.

A Historical-Genetic Reconstruction of Human Extra-Pair Paternity.

Paternity testing using genetic markers has shown that extra-pair paternity (EPP) is common in many pair-bonded species [1, 2]. Evolutionary theory and empirical data show that extra-pair copulations can increase the fitness of males as well as females [3, 4]. This can carry a significant fitness cost for the social father, who then invests in rearing offspring that biologically are not his own [5]. In human populations, the incidence and correlates of extra-pair paternity remain highly contentious [2, 6, 7]. Here, we use a population-level genetic genealogy approach [6, 8] to reconstruct spatiotemporal patterns in human EPP rates. Using patrilineal genealogies from the Low Countries spanning a period of over 500 years and Y chromosome genotyping of living descendants, our analysis reveals that historical EPP rates, while low overall, were strongly impacted by socioeconomic and demographic factors. Specifically, we observe that estimated EPP rates among married couples varied by more than an order of magnitude, from 0.4% to 5.9%, and peaked among families with a low socioeconomic background living in densely populated cities of the late 19th century. Our results support theoretical predictions that social context can strongly affect the outcomes of sexual conflict in human populations by modulating the incentives and opportunities for engaging in extra-pair relationships [9-11]. These findings show how contemporary genetic data combined with in-depth genealogies open up a new window on the sexual behavior of our ancestors.

Estimating Y-Str Mutation Rates and Tmrca Through Deep-Rooting Italian Pedigrees.

In the population genomics era, the study of Y-chromosome variability is still of the greatest interest for several fields ranging from molecular anthropology to forensics and genetic genealogy. In particular, mutation rates of Y-chromosomal Short Tandem Repeats markers (Y-STRs) are key parameters for different interdisciplinary applications. Among them, testing the patrilineal relatedness between individuals and calculating their Time of Most Recent Common Ancestors (TMRCAs) are of the utmost importance. To provide new valuable estimates and to address these issues, we typed 47 Y-STRs (comprising Yfiler, PowerPlex23 and YfilerPlus loci, the recently defined Rapidly Mutating [RM] panel and 11 additional markers often used in genetic genealogical applications) in 135 individuals belonging to 66 deep-rooting paternal genealogies from Northern Italy. Our results confirmed that the genealogy approach is an effective way to obtain reliable Y-STR mutation rate estimates even with a limited number of samples. Moreover, they showed that the impact of multi-step mutations and backmutations is negligible within the temporal scale usually adopted by forensic and genetic genealogy analyses. We then detected a significant association between the number of mutations within genealogies and observed TMRCAs. Therefore, we compared observed and expected TMRCAs by implementing a Bayesian procedure originally designed by Walsh (2001) and showed that the method yields a good performance (up to 96.72%), especially when using the Infinite Alleles Model (IAM).

A game of hide and seq: Identification of parallel Y-STR evolution in deep-rooting pedigrees.

Short tandem repeats on the Y-chromosome (Y-STRs) are common DNA polymorphisms useful for genetic genealogy, population and evolutionary genetics, human genetics, pathology and forensic sciences. It is important to identify all Y-STR variants and to have knowledge of Y-STR mutation rates in order to correctly estimate the time to the most recent common ancestor (tMRCA) between paternally related individuals. When capillary electrophoresis (CE) is performed to analyze genealogical pairs, Y-STR sequence variations remain hidden when the number of repeats is identical. These hidden variations could be due to parallel Y-STR changes or modifications (PM) that occur independently in different lineages leading to alleles with identical number of repeats. In this study, we detect for the first time twelve PM by analyzing 133 males (960 meiosis) in extended deep-rooting family pedigrees on 42 Y-STRs. These PM were observed in nine Y-STR loci with mutation rates of at least 5.94 × 10-3 per generation. Sequencing analysis made it possible to distinguish insertions/deletions in different repeat regions revealing the presence of two unique changes in three PM on rapidly mutating and complex Y-STRs DYS724-ab and DYS518. Sequencing unraveled more information concerning the identity of alleles, and increased allelic discrimination possibilities which is of great importance in population genetics and forensic analysis. Limiting the analysis to CE could lead to wrong ancestral allele assumptions, to false negative interpretations and to tMRCA underestimations. These observations highlight the importance and added value of sequencing analysis and suggest a shift in genotyping methods from CE to next generation sequencing.

Spatially Dense 3D Facial Heritability and Modules of Co-heritability in a Father-Offspring Design.

Introduction: The human face is a complex trait displaying a strong genetic component as illustrated by various studies on facial heritability. Most of these start from sparse descriptions of facial shape using a limited set of landmarks. Subsequently, facial features are preselected as univariate measurements or principal components and the heritability is estimated for each of these features separately. However, none of these studies investigated multivariate facial features, nor the co-heritability between different facial features. Here we report a spatially dense multivariate analysis of facial heritability and co-heritability starting from data from fathers and their children available within ALSPAC. Additionally, we provide an elaborate overview of related craniofacial heritability studies. Methods: In total, 3D facial images of 762 father-offspring pairs were retained after quality control. An anthropometric mask was applied to these images to establish spatially dense quasi-landmark configurations. Partial least squares regression was performed and the (co-)heritability for all quasi-landmarks (∼7160) was computed as twice the regression coefficient. Subsequently, these were used as input to a hierarchical facial segmentation, resulting in the definition of facial modules that are internally integrated through the biological mechanisms of inheritance. Finally, multivariate heritability estimates were obtained for each of the resulting modules. Results: Nearly all modular estimates reached statistical significance under 1,000,000 permutations and after multiple testing correction (p ≤ 1.3889 × 10-3), displaying low to high heritability scores. Particular facial areas showing the greatest heritability were similar for both sons and daughters. However, higher estimates were obtained in the former. These areas included the global face, upper facial part (encompassing the nasion, zygomas and forehead) and nose, with values reaching 82% in boys and 72% in girls. The lower parts of the face only showed low to moderate levels of heritability. Conclusion: In this work, we refrain from reducing facial variation to a series of individual measurements and analyze the heritability and co-heritability from spatially dense landmark configurations at multiple levels of organization. Finally, a multivariate estimation of heritability for global-to-local facial segments is reported. Knowledge of the genetic determination of facial shape is useful in the identification of genetic variants that underlie normal-range facial variation.

Recombination hotspots in an extended human pseudoautosomal domain predicted from double-strand break maps and characterized by sperm-based crossover analysis.

The human X and Y chromosomes are heteromorphic but share a region of homology at the tips of their short arms, pseudoautosomal region 1 (PAR1), that supports obligate crossover in male meiosis. Although the boundary between pseudoautosomal and sex-specific DNA has traditionally been regarded as conserved among primates, it was recently discovered that the boundary position varies among human males, due to a translocation of ~110 kb from the X to the Y chromosome that creates an extended PAR1 (ePAR). This event has occurred at least twice in human evolution. So far, only limited evidence has been presented to suggest this extension is recombinationally active. Here, we sought direct proof by examining thousands of gametes from each of two ePAR-carrying men, for two subregions chosen on the basis of previously published male X-chromosomal meiotic double-strand break (DSB) maps. Crossover activity comparable to that seen at autosomal hotspots was observed between the X and the ePAR borne on the Y chromosome both at a distal and a proximal site within the 110-kb extension. Other hallmarks of classic recombination hotspots included evidence of transmission distortion and GC-biased gene conversion. We observed good correspondence between the male DSB clusters and historical recombination activity of this region in the X chromosomes of females, as ascertained from linkage disequilibrium analysis; this suggests that this region is similarly primed for crossover in both male and female germlines, although sex-specific differences may also exist. Extensive resequencing and inference of ePAR haplotypes, placed in the framework of the Y phylogeny as ascertained by both Y microsatellites and single nucleotide polymorphisms, allowed us to estimate a minimum rate of crossover over the entire ePAR region of 6-fold greater than genome average, comparable with pedigree estimates of PAR1 activity generally. We conclude ePAR very likely contributes to the critical crossover function of PAR1.

The biological relevance of a medieval king's DNA.

The discovery of the presumably lost grave of the controversial English king Richard III in Leicester (U.K.) was one of the most important archaeological achievements of the last decennium. The skeleton was identified beyond reasonable doubt, mainly by the match of mitochondrial DNA to that of living maternal relatives, along with the specific archaeological context. Since the genetic genealogical analysis only involved the DNA sequences of a single 15th century individual and a few reference persons, biologists might consider this investigation a mere curiosity. This mini-review shows that the unique context of a historical king's DNA also has relevance for biological research per se - in addition to the more obvious historical, societal and educational value. In the first place, the historical identification appeared to be a renewed forensic case realising a conservative statement with statistical power based on genetic and non-genetic data, including discordant elements. Secondly, the observation of historical non-paternity events within Richard III's patrilineage has given rise to new research questions about potential factors influencing the extra-pair paternity rate in humans and the importance of biological relatedness for the legal recognition of a child in the past. Thirdly, the identification of a named and dated skeleton with the known historical context serves as a reference for bioarchaeological investigations and studies on the spatio-temporal distribution of particular genetic variance. Finally, the Richard III case revealed privacy issues for living relatives which appear to be inherent to any publication of genetic genealogical data.

Mediterranean Y-chromosome 2.0-why the Y in the Mediterranean is still relevant in the postgenomic era.

CONTEXT: Due to its unique paternal inheritance, the Y-chromosome has been a highly popular marker among population geneticists for over two decades. Recently, the advent of cost-effective genome-wide methods has unlocked information-rich autosomal genomic data, paving the way to the postgenomic era. This seems to have announced the decreasing popularity of investigating Y-chromosome variation, which provides only the paternal perspective of human ancestries and is strongly influenced by genetic drift and social behaviour. OBJECTIVE: For this special issue on population genetics of the Mediterranean, the aim was to demonstrate that the Y-chromosome still provides important insights in the postgenomic era and in a time when ancient genomes are becoming exponentially available. METHODS: A systematic literature search on Y-chromosomal studies in the Mediterranean was performed. RESULTS: Several applications of Y-chromosomal analysis with future opportunities are formulated and illustrated with studies on Mediterranean populations. CONCLUSIONS: There will be no reduced interest in Y-chromosomal studies going from reconstruction of male-specific demographic events to ancient DNA applications, surname history and population-wide estimations of extra-pair paternity rates. Moreover, more initiatives are required to collect population genetic data of Y-chromosomal markers for forensic research, and to include Y-chromosomal data in GWAS investigations and studies on male infertility.